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Abstract 

Thermal fluctuations of the granular gas under the homogeneous cooling state 
(HCS) are estimated using two-point kinetic theory by Tsuge-Sagara. Thermal fluc- 
tuations of the elastic gas are modified for the granular gas by nonequilibrium mo- 
ments, which defines the distribution function under the HCS. The deviations of 
thermal fluctuations for the granular gas from those for the elastic gas obtained by 
the fluctuation-dissipation theorem are calculated as a function of the restitution 
coefficient. 

PACS numbers: 



I. INTRODUCTION 



Q. 



The granular gas becomes a popular issue owing to its rich physical phenomena [1[ . Especially, the 
pattern formation via the inelastic collisions is presumably caused by the local thermal fluctuations 
and fluid dynamic fluctuations. The decay of the amplitude of fluctuations via the decrease of the 
temperature as a result of inelastic collisions yields a larger pattern formation via the dissipation. On 
the other hand, properties of thermal fluctuations for the granular gas are not understood completely, 
whereas Brey et al. [2| formulated the thermal fluctuations on the basis of the the linear Boltzmann 
equation with a fluctuating white-noise term. In this paper, we consider on deviations of thermal 
fluctuations for the granular gas from those for the elastic gas. To simplify our discussions, we restrict 



a. 



ourselves to the HCS [2jj. To evaluate thermal fluctuations for the granular gas under the HCS, 
the two-point kinetic theory by Tsuge-Sagara [s| for the elastic gas is extended into the granular gas. 
Assuming the HCS, the correlation between different particles, which yields fluid dynamic fluctuations 
{jj, can be neglected, and only self-correlation function [jj] is considered exclusively. Finally, deviations 
of thermal fluctuations for the granular gas from those for the elastic gas are obtained as a function 
of nonequilibrium moments, which define the distribution function under the HCS. In this letter, the 
two-point kinetic theory is revisited in detail, and extended into the granular gas. Finally, we evaluate 
deviations of thermal fluctuations for the granular gas from those for the elastic gas as a function of 
the restitution coefficient. In this paper, inelastic hard sphere particles with a constant restitution 
coefficient and the smooth surface are discussed. 



II. REVISITED TWO-POINT KINETIC THEORY BY TSUGE-SAGARA 



,3] 



The two-point kinetic theory is proposed by Tsuge-Sagara, and discussed for the elastic gas. In this 
section, we revisit the two-point kinetic theory and introduce fluctuation-dissipation (FD) theorem 4] 
for thermal fluctuations of the elastic gas. The stochastic Boltzmann equation is the most fundamental 
equation in the two-point kinetic theory. The stochastic Boltzmann equation is written as follows: 



0, (1) 



ac) [f?(t,c)f?(t,c)" 



where g(t, Q is the microscopic distribution function defined by: 

N 

e(t,C) = £ tf [c-C w (*)]> (3) 

s=l 

where s is the index of the partiles, N is the total number of particles, £ = (c, x), in which c E {Y 3 : 
— oc < c < oo} is the velocity space and x is the physical space. In eq. (1), the prime indicates the 
state after binary collisions, g = \c — c\ is the relative velocity of the binary colliding particles, a is 
the differential collision-cross section, x is the deflection angle, e is the scattering angle. 
For convenience, we set £ = (t, £). In £-space, we consider the correlation between two points "a" and 
"/3". As a result, g(a) = g(£(a)) and g{(3) = g(£((3)) are used. At point a, the stochastic Boltzmann 



3 



equation in eq. (1) is rewitten as: 

d 



B(a)g (a) 



dt(a] 



+ c (a, 







dx(a) 



g(a) — J* \6l\ol) \g{a)g{a)\ = 0. 



(4) 



Next, we define the averaged quantites in £-space as follows using finite volumes AV G x and AV C 6 



1 



dx / dap, 



^ AV C AV 

Aip = ip — ip, 



(5) 
(6) 



where ip is the arbitrary function defined in £-space. 

On the basis of eqs. (5) and (6), we obtain the following equation by multiplying Ag(j3) by both sides 
of eq. (4) and taking average in £-space, 



d . . d 

+ c(a) 



dt(a) 



dx(a) 



J (a\a) g(a)g(a)g(P) - g(a) g(a) f ((3) 



(7) 



where f(a) = g(a) A f{(3) = g((3). 



In eq. (7), g(a)g(f3) is decomposed as: 



g(a)g(/3) = f n {a; f3) + g(a; f5), 



(8) 



where /// (a; /3) is the two-point phase density of different particles and g(a; (3) is the probability of 
finding same particles at t = t(a) A £ = C( a ) an d t = t((3) A £ = C(P)- fn{ a ]P) is related to the 
hydrodynamic fluctuation term (p(a; (5) by 



fn(a;P) = + (1 + 



(a; /3) +/(«)/(/?), (For 1 « iV), 



(9) 
(10) 



when 1 CJV and i(a) — >■ t((3), we obtain 



lim Jn(a;P) = <Ka,l3) + f(a)f(P), 
lim o(a;0)=*[£(a)-£O3)]/(a) 



(11) 
(12) 



Assuming that the hydrodynamic fluctuations is ignorable, (i.e., 4>(ct; f3) = or fn(a; j3) = f(a)f(/3)), 
we obtain the following equation from eq. (7): 

d 



g (a; (3) = / (a\a) (f(a)g(a; (3) + f(a)g(a; £)) , 



(13) 
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where r = t(a) — t(B). 

n 

Tsuge-Sagara expanded g(a; j3) using the Grad's method [5j as follows: 



g(a;P) = u)(a)u)(/3) ^ 



c J + K J\K\ l i 



(14) 



where c = V i?T, in which is the gas constant and T is the temperature, oj(a) = j^c 1 ) 2 ex P 



and H^p is the Hermite polynomial 



c(q) 2 



^ ij...lm.. 



(((a) - C(/3),r) = | c J+K H\p {a)H^ (f3)g(a, /3)dv(a)dv(f3). (15) 
We assume the initial equilibrium condition for g(a,/3) using eq. (12) as follows 



b(«,/3)] T=0 = 5[e(«)-e(/3)]/MB(a), 

where /mb is the Maxwell-Boltzmann distribution function. 
Consequently, Q is; 

QiSm.M") ~ = nc J+K 5[i(a) - £(/3)]<5 JK V 



(16) 



(17) 



In particular, we have the correlation of the fluctuation of the number density and the correlation of 
the fluctuation of the velocity at r = as follows: 



[£(0,0)- 



T=0 



T=0 



nS [x(a) — x((3)] 

(correlation of fluctuation of number density), 

n~ 1 c 2 5u5 [x(a) — x((3)] 

(correlation of fluctuation of velocity) . 



(18) 



(19) 



In a similar way, we have the correlation of the fluctuation of the static pressure, the correlation of 
the fluctuation of the pressure deviator, and the correlation of the fluctuation of the heat flux at r = 
as follows: 



Ap(a)Ap(/3) 



,(2,2) 
Hj,lm 



(3,3) 
l i,3 



T = 



T = 



T = 



p 2 n~ 1 (5/3)5[x(a) - x(P)} 

(correlation of fluctuation of static pressure) 

P 2 ( 2 

—5 [x(a) - x{(3)\ ( 5u5 jm + 5ji5 im - -5ij5i m 

ft \ O 

(correlation of fluctuation of pressure deviator) 
2~n 



RT5 [x(a) - x(p)] 6^ 



(correlation of fluctuation of heat flux), 



(20) 



(21) 



(22) 



where v = nkT is the static pressure. 



13 moment equations 



are obtained by neglecting the convective terms from eq. (13) as follows: 



#Q (0 ' 0) = ^Qfi l) = ^Q (2 < 2) = o, 

OT OT ' OT 

9 (2,2) „ „ (2,2) „ 
Q-%\lL + 6nS ?^Zm = °> 

9 (3,3) . . „ (3,3) n 

111 + 4nB <; = 0, 



where 6nB = p/n, in which [i is the viscosity coefficient for the elastic gas. 
As a result, we obtain; 

ifjM = n~ 1 p 2 exp (-rp/fj,) 5 [x(a) - x((3)] ^5u5 jm + 6 im 8ji - ^8^5^ 
qf& = n-VC p Texp {-rpCp/X) [x(a) - x{p)] 5 iU 



(23) 
(24) 
(25) 



(26) 
(27) 



where C p = 5R/2 is the isobaric specific heat and A is the thermal conductivity of the elastic gas. 
By approximating e~ aT ~ (2/a)5(r) for 1 <^ a, we obtain thermal fluctuations obtained by the FD 
theorem 4 ] . 



III. THERMAL FLUCTUATIONS OF GRANULAR GAS USING TWO-POINT 

KINETIC THEORY 

The two-point kinetic theory of the elastic gas by Tsuge-Sagara was revisited. On the basis of 
forgoing discussions, we consider thermal fluctuations for the granular gas. 
Different points from the elastic gas are; 

1. The collision term defined in eqs. (1) and (2) is replaced by the inelastic Boltzmann equation. 

2. The distribution function of the granular gas under the HCS is not the Maxwell-Boltzmann 
distribution. 

Consequently, the initial self-correlation function [g(a, /3)] T=Q = <5(£(a) — C(/5))/ ( - ^ is not ex- 
pressed by the Maxwell-Boltzmann distribution as shown in eq. (16). 

To simplify our discussion, the restitution coefficient of the granular gas is constant in regardless of 
the relative velocity g. From the item "1", the collision term of the stochastic Boltzmann equation in 



eq. (2) is replaced by 



Q: 



r n (cic) [ e <t,Q e (t,t) 



27T 



1 



(28) 



Q. 



where the double prime indicates the inverse collision e is the constant restitution coefficient and 
superscript in indicates the inelastic collisions of the granular gas. 

Under the HCS, we can start our discussion from eq. (13) by replacing the collision term J! by 
in eq. (28). As discussed in the item "2", the distribution function /(°) under the HCS does not obey 
the Maxwell-Boltzmann distribution function. Generally, /(°) under the HCS is expanded around the 
Maxwell-Boltzmann distribution as follows: 



/(°) = f MB ( 1 + l - a ^H^ + -±-a^H^ + ■■■ + 9 2k a 2k H^ + • 
J V 8 5040 



(29) 



where a^ 2k ^ is the spherically symmetric Grad's moment 8] and #( 2fc ) is the coefficient for a normal- 
ization. 

For example, and are written as follows 



ff( 4 ) = v ± - l 0u 2 + 15, 

H (6) = v 6 - 2lv A + 105v 2 - 105, 



(30) 
(31) 



where v 2 = c 2 /(RT). 

On the basis of eq. (29), eqs (18)-(22) can be extended into the case of granular gas as follows: 



[Q m }% 



T=0 

in 



r=0 



(2j2)1 m 



(3,3) 



T=0 

in 

T=0 



n5 [x(a) — x(f3)\ , 
n- 1 c 2 5 il 5[x(a)-x(f3)}, 

x (5/3) (lW 4) )<5 [x(a) — x(/3)] , 
^-<5 [x(a) - x(/3)] (l + a (4) ) (s a 8j m + 5jiS im - ^SijSiJ] 

1 + ^ a (4) + l a (fi)\ P^ RT5 [x{a) _ x{p)] S 

2 4 12 J n 1 y ' y n 



(32) 
(33) 
(34) 
(35) 
(36) 
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Moment equations for the elastic gas in eqs. (23)- (25) are extended into the granular gas as follows: 



d 
dr 


q(0,0) 


in Q 

= 




in 


= 0, 




d 
dr 




in 

= -c 


Q(2,2)" 


in 






d 
dr 


' (2,2)" 
Qij,lm 




D 


" (2,2) 
^ij,lm 


in 

= 


d 


- (3,3)" 
%l 


m + 4n^ (3 ' 3) (e) 


B 


_(3,3)- 
^i,l 


in 

= 0, 



(37) 
(38) 
(39) 
(40) 



where £ is the cooling rate 

^( 2 ' 2 )(e) and J^( 3 ' 3 ) (e) in eqs. (39) and (40) are given by Jenkins- Richman 6] as follows: 



^(2,2) (t 



(l + 6)(3-e) 



^(3,3) (e) = (l_+e)(49-33e) 



32 



(41) 
(42) 



As a result of eqs. (35), (36), (39) and (40), we obtain 

1 + a 



■ (2,2) l m 
Qij,lm 



ew{-Ty {2 ' 2) (e)p/fi) (2i2) 



exp(-rp/^i) 



,{3,3)1 1 

hi 



1 + ii a (4) + l a (6) 

2 30 



exp(-rPr^ 3 )(6W/i) (3j3) 

1i ,l ' 



(43) 
(44) 



exp (— r Prp//i) 
where Pr = 2/3 is performed for the monatomic gas. 

By approximating e~ ar ~ (2/a)<5(r) for 1 <C a, we obtain the deviations of thermal fluctuations of the 
granular gas from those for the elastic gas obtained by the FD theorem as follows: 

1 



(2,2) y n 



(3,3) 



^(2,2) ( £ ) 



1 + a { 



m „(2,2) 

lij,lm' 



5(2,2) 



111 



1 



1 + 11 (4) + Jl q (6) \ 0,3) 

^(3,3) ( e ) V 2 a + 30° 

5(3,3) 



(45) 



(46) 



where S , ( 2,2 ) is the function, which reveals the deviation of thermal fluctuations of the pressure deviator 
for the granular gas from those for the elastic gas, and S^ 3,3 -* is the function, which reveals the deviation 
of thermal fluctuations of the heat flux for the granular gas from those for the elastic gas. 



IV. DISCUSSION 



Finally, and in eqs. (45) and (46) are considered in detail, and under the HCS 

is calculated by Brilliantov-Poschel as follows 7(: 



«( 4 ) = 

-16(-1623 + 1934e + 895e 2 - 364e 3 + 3510e 4 - 7424e 5 + 3321e 6 - 480e 7 + 240e 8 ) 
214357 - 172458e + 112155e 2 + 25716e 3 - 4410e 4 - 84480e 5 + 34800e 6 - 5600e 7 + 2800e 8 ' 

13440(217 - 386e - 669e 2 + 1548e 3 + 154e 4 - 1600e 5 + 816e 6 - 160e 7 + 80e 8 ) 
214357 - 172458e + U2155e 2 + 25716e 3 - 4410e 4 - 84480e 5 + 34800e 6 - 5600e 7 + 2800e 8 ' 



(47) 



From eqs. (41), (42), (45), (46) and (47), 5^ 2 ' 2) and S {3 ^ are plotted in Fig. 1. 



S (2 ' 2) and S (3 ' 3) versus £ 




(2,2) 
(3,3) 



FIG. 1 and versus e (restitution cofficient) 

5 (2,2) < 1 at 0.786 < e < 1 and 1 < 5 (2 ' 2) at < e < 0.786 as shown in Fig. 1. 5 (3,3) < 1 at 
0.36 < e < 1 and 1 < 5 (3,3) at < e < 0.36 as shown in Fig. 1. S"( 2 ' 2 ) has its minimum value 
= 0.986 at e = 0.894 and has its minimum value sQ® = 0.731 at e = 0.749. The 

deviation- function jS^ 2,2 ) does not change from unity around e ~ 1. As a result, the formulation of 
thermal fluctuations of the pressure deviator for the granular gas is obtained by slightly modifying that 
for the elastic gas in weakly inelastic regime (i.e., e ~ 1). On the other hand, the deviation-function 
,5(3,3) drastically changes from unity around e ~ 1. As a result, the formulation of thermal fluctuations 
of the heat flux for the granular gas is obatined by remarakably modifying that for the elastic gas even 



9 



in the weakly inelastic regime. 
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